The Role of Gorse (Ulex parviflorus Pourr. Scrubs) in a Mediterranean Shrubland Undergoing Climate Change: Approach by Hyperspectral Measurements

The goal of this study was to observe the neighbor effect of Gorse, a plant of the Fabaceae family, on three typical species of Mediterranean shrubland: kermes oak, white Cistus and rosemary. For this purpose, a hyperspectral analysis and the application of vegetation indices (VIs) were carried out. We provide the spectral signature of Gorse, which differs mainly from that of its companion species in the band between 700 and 1350 nm. This supposed Gorse effect was tested in natural conditions and in conditions of forced drought to simulate the effects of the climate change predicted for the Mediterranean Basin. Field spectrometry demonstrated the existence of such interactions between the four species. In control stands, the presence of Gorse significantly modifies the spectral responses of kermes, white Cistus and rosemary, mainly in the near-infrared region (700–1350 nm). Both tri- and tetra-specific plant assemblages also exhibited spectral changes, suggesting an indirect effect of Gorse. Under drought conditions, one-way ANOVA followed by Fisher’s LSD test led us to identify the features involved in plants’ coexistence with Gorse. The Cistus albidus reflectance spectrum was clearly increased in the presence of Gorse in rain-exclusion conditions. The application of several VIs allowed us to extract new information on the variation of spectral signatures. Unexpectedly, nitrogen supply by Gorse was not shown, except for Cistus, as shown by the VI NDVI (N) analysis. However, this study proved that Gorse can modify the behavior of its companion species in controls, but also in drought conditions, by increasing their photosynthesis activity (NIRvP) and water content (ratio R975/R900). Gorse therefore appears as a key species in the ecosystem of the Mediterranean shrubland, but its high vulnerability to drought leaves a vacant ecological niche in plant communities. While the spectral reflectance increases linearly with the specific richness in the lack of any disturbance, by contrast, climate aridification imposes a double reciprocal profile. This clearly means that multispecific plant communities cope better with climate change. Nevertheless, knowledge of the underlying mechanisms requires further structural, chemical, and biochemical investigation.


Introduction
Understanding the processes that govern plant-community structure and composition is one of the major challenges in ecology. Ulex parviflorus Pourr. Scrubs (Gorse) belongs to the family Fabaceae, which is widespread in Mediterranean shrublands. As a legume, improving soil-nutritional quality [1,2], it is also described as a fire-prone species, in accordance with the 'kill thy neighbour' hypothesis [3,4]. These plant behaviors, in terms of both atmospheric-nitrogen fixation and flammability, are important drivers in terrestrial ecosystems and, thus, play a fundamental role in species evolution and ecosystem dynamics. Provence). Two types of device are implanted in plots of 16 m². Rain-exclusion devices (Figure 1A) equipped with gutters and control devices with inverted gutters are used to consider the shading effect ( Figure 1B). These devices are installed on plots with different natural combinations of species; a total of 93 plots are arranged on the CLIMED site. For this study, several types of plant combination were considered (11 assemblies). Monospecific plots for the 4 dominant woody species: Ulex parviflorus (Up), Cistus albidus (Ca), Quercus coccifera (Qc), Rosmarinus officinalis (Ro). Bispecific plots with the following associations: UpCa, UpQc, UpRo. Trispecific plots with the following plant combinations: QcCaUp, QcRoUp, CaRoUp, CaQcUp, RoQcUp, RoCaUp. Tetraspecific plots: QcCaRoUp, CaQcRoUp and RoQcCaUp. The first name of species of each assemblage corresponded to the species that was considered for spectral measurements. Each modality of water treatment (control and exclusion devices) was repeated 3 times. Thus, a total of 66 plots were selected to test the neighbor effects of U. parviflorus.

Selection of Individuals
Plants were chosen central to the plot to avoid border effect. The selected individuals were located at least 1 m from the edge of the plot. In each plot, 3 individuals were carefully selected; they were located as close as possible to the other species for maximum interaction.

Spectral Measurements
The measurements were performed during a field campaign (27 April 2022 to 10 May 2022). The reflectance was obtained between 350 and 2500 nm using a portable spectroradiometer (FieldSpec ® 4 ASD, USA) with a field of view of 25°. First, a reflectance calibration was performed using a Spectralon ® disk before plant measurement. For the measurements, the optical fiber was positioned at nadir, 10 cm above the canopy of target plant. On each of the 3 individuals per plot, 3 reflectance measurements were performed. To check for sensor stability, each measurement was repeated 3 times. The data acquisition was performed between 11:00 and 15:00 in similar synoptic conditions. A total of 1215 data were recorded: monospecific (189); bispecific (324); trispecific (486); and tetraspecific (216).

Spectral Post-Treatment
With ViewSpecPro, spectra were corrected by applying a "splice correction" for linear interpolation at the spectral discontinuities. On all data obtained, artifacts representing absorption bands of atmospheric water were observed. Therefore, the bands from 1340 nm to 1430, from 1810 nm to 1980 nm and above 2400 nm were removed for all the data. For this study, several types of plant combination were considered (11 assemblies). Monospecific plots for the 4 dominant woody species: Ulex parviflorus (Up), Cistus albidus (Ca), Quercus coccifera (Qc), Rosmarinus officinalis (Ro). Bispecific plots with the following associations: UpCa, UpQc, UpRo. Trispecific plots with the following plant combinations: QcCaUp, QcRoUp, CaRoUp, CaQcUp, RoQcUp, RoCaUp. Tetraspecific plots: QcCaRoUp, CaQcRoUp and RoQcCaUp. The first name of species of each assemblage corresponded to the species that was considered for spectral measurements. Each modality of water treatment (control and exclusion devices) was repeated 3 times. Thus, a total of 66 plots were selected to test the neighbor effects of U. parviflorus.

Selection of Individuals
Plants were chosen central to the plot to avoid border effect. The selected individuals were located at least 1 m from the edge of the plot. In each plot, 3 individuals were carefully selected; they were located as close as possible to the other species for maximum interaction.

Spectral Measurements
The measurements were performed during a field campaign (27 April 2022 to 10 May 2022). The reflectance was obtained between 350 and 2500 nm using a portable spectroradiometer (FieldSpec ® 4 ASD, USA) with a field of view of 25 • . First, a reflectance calibration was performed using a Spectralon ® disk before plant measurement. For the measurements, the optical fiber was positioned at nadir, 10 cm above the canopy of target plant. On each of the 3 individuals per plot, 3 reflectance measurements were performed. To check for sensor stability, each measurement was repeated 3 times. The data acquisition was performed between 11:00 and 15:00 in similar synoptic conditions. A total of 1215 data were recorded: monospecific (189); bispecific (324); trispecific (486); and tetraspecific (216).

Spectral Post-Treatment
With ViewSpecPro, spectra were corrected by applying a "splice correction" for linear interpolation at the spectral discontinuities. On all data obtained, artifacts representing absorption bands of atmospheric water were observed. Therefore, the bands from 1340 nm to 1430, from 1810 nm to 1980 nm and above 2400 nm were removed for all the data.

Spectral Indices
Vegetation Indices (VIs) are widely used to predict biophysical properties of natural ecosystems. Their principal advantage is that they enhance spectral-reflectance information by detecting the spectral variability that might occur, for instance, in plant canopies or physiological and structural properties of leaves [16]. First, we considered two VIs that are related to plant-water content and photosynthesis.
Water-content index was calculated as follows [17]. R975/R900. NIRvP: structural proxy for SIF (solar-induced fluorescence). Photosynthesis [18]: NIRvP = NIRv × PAR, where PAR is the photosynthetically active radiation. The NIR is the average reflectance between 770 and 780 nm. NDVI = N IR−RED re f N IR+RED re f , where RED ref is the average reflectance between 650 and 660 nm. The PAR were collected from the meteorological station of the CLIMED site at each time of canopy-reflectance measurement. Second, other VIs were selected to assess plant nitrogen content, as well as their structure and functioning either directly or indirectly linked to photosynthesis. These are recorded in Table 1.

Statistical Analysis
Spectral data were processed from the Metaboanalyst platform [31]. They were first centered and reduced, after which they were normalized for multivariate analysis (partial least-square discriminant analysis (PLSDA)), followed by ANOVA and Fisher's LSD post hoc tests. The latter allows pairwise comparison of assemblages for each wavelength of spectrum. In the same way, spectral indices were processed by PLSDA and then classified by the VIP scores (importance of the variables in projection on axis 1 of PLSDA). The most discriminatory indices that fitted the objectives of the study were processed by twoway-ANOVA with Statgraphics. Correlation models between total reflectance and species richness were also created with Statgraphics.

Scanning Electronic Microscopy (SEM) Observations
Three twigs of Ulex parviflorus, (3 individuals) at different phenological stages (young and mature) were observed with a Zeiss scanning electron microscope, SEM EVO 15 (with lanthanum boride source), in wet mode.

Spectral Signature of Ulex Parviflorus
The PLSDA analysis exhibited a total inertia of 37% in accordance with a weak discrimination between Gorse, kermes oak and rosemary ( Figure 2). However, U. parviflorus seems spectrally closer to kermes than to rosemary and cistus. Exploring the spectral profiles showed differences in some areas of the spectrum (Figure 3). The Gorse sems to have exhibited much lower reflectance than the other species in the band between 700 and 1400 nm. This was confirmed by the one-way-ANOVA, followed by the Fisher's LSD test, with the identification of the wavelengths that significantly distinguished the spectrum of U. parviflorus from the 3 other species (Supplementary Table S1).
ANOVA with Statgraphics. Correlation models between total reflectance and species richness were also created with Statgraphics.

Scanning Electronic Microscopy (SEM) Observations
Three twigs of Ulex parviflorus, (3 individuals) at different phenological stages (young and mature) were observed with a Zeiss scanning electron microscope, SEM EVO 15 (with lanthanum boride source), in wet mode.

Spectral Signature of Ulex Parviflorus
The PLSDA analysis exhibited a total inertia of 37% in accordance with a weak discrimination between Gorse, kermes oak and rosemary ( Figure 2). However, U. parviflorus seems spectrally closer to kermes than to rosemary and cistus. Exploring the spectral profiles showed differences in some areas of the spectrum (Figure 3). The Gorse sems to have exhibited much lower reflectance than the other species in the band between 700 and 1400 nm. This was confirmed by the one-way-ANOVA, followed by the Fisher's LSD test, with the identification of the wavelengths that significantly distinguished the spectrum of U. parviflorus from the 3 other species (Supplementary Table S1).

Hyperspectral Response of Interspecific Assemblages
The spectral data obtained according to the different types of assemblage tested under natural conditions (control modality) showed that the bands from 700 to 1400 nm were

Hyperspectral Response of Interspecific Assemblages
The spectral data obtained according to the different types of assemblage tested under natural conditions (control modality) showed that the bands from 700 to 1400 nm were the most discriminatory ( Figure 4). The direct effect of Gorse on its companion species may be observed by comparing the spectra obtained in monospecific and bispecific stands. The pairwise comparison of the assemblages for each wavelength resulted in more than one thousand significant differences (p < 0.05). We present in this report the most significant differences obtained for each species ( Figure 5). The Gorse did not directly affect the reflectance of Ca at 1030 nm (Figure 5a), nor that of Ro at 559 nm ( Figure 5b). However, an increase in reflectance occurred for QcUp compared to Qc at 1291 nm ( Figure 5c). Conversely, all the three species (UpCa; UpRo; UpQc) contributed to the increase in the reflectance of Gorse at 2254 nm compared to their monospecific stand ( Figure 5d). By plotting the integrated reflectance with the specific richness, we found a positive and linear correlation, meaning that the spectral reflectance increased with species richness (Supplementary Figure S1).    Figure 6 compares the reflectance of each assemblage as a function of the control and exclusion treatment modalities. The features identified by the one-way ANOVA followed by the Fisher's LSD test are summarized in Supplementary Table S2. In general, there was no significant effect of rainfall exclusion on the monospecific stands (Figure 6a,e,i, Supplementary Table S2). In the bispecific stands, there was no significant effect of Up on Ro (RoUp); this effect was only observed in the SWIR2 band (1801-2350 nm). No difference was found for QcUp all over the spectrum (Supplementary Table S2). The most significant effects were observed between eCaUp and cCaUp (Supplementary Table S2), where the reflectance was 12% higher than the cCaUp at its maximum value, about 1140 nm ( Figure  6f). Up appears to play a role in the reflectance capacity of Ca under drought conditions.  Figure 6 compares the reflectance of each assemblage as a function of the control and exclusion treatment modalities. The features identified by the one-way ANOVA followed by the Fisher's LSD test are summarized in Supplementary Table S2. In general, there was no significant effect of rainfall exclusion on the monospecific stands (Figure 6a,e,i, Supplementary Table S2). In the bispecific stands, there was no significant effect of Up on Ro (RoUp); this effect was only observed in the SWIR2 band (1801-2350 nm). No difference was found for QcUp all over the spectrum (Supplementary Table S2). The most significant effects were observed between eCaUp and cCaUp (Supplementary Table S2), where the reflectance was 12% higher than the cCaUp at its maximum value, about 1140 nm (Figure 6f). Up appears to play a role in the reflectance capacity of Ca under drought conditions. For the trispecific combination, the lowest reflectance was recorded in the  (Figure 6g,k). This decrease in reflectance was also observed in the tetraspecific eRoQcCaUP and eQcCaRoUp stands. For the effect of the three species on Up, we noted an inverse response for the tetraspecific assemblage (Figure 6p). In general, the reflectance in the exclusion decreased with the species richness. This property did not follow a linear pattern, but a double inverse profile (Supplementary Figure S2). For the trispecific combination, the lowest reflectance was recorded in the exclusion, especially from the eCaQcUp and eRoQcUP stands (Figure 6g,k). This decrease in reflectance was also observed in the tetraspecific eRoQcCaUP and eQcCaRoUp stands. For the effect of the three species on Up, we noted an inverse response for the tetraspecific assemblage (Figure 6p). In general, the reflectance in the exclusion decreased with the species richness. This property did not follow a linear pattern, but a double inverse profile (Supplementary Figure S2).

Water-Content Index: The Ratio R975/R900
High reflectance in water absorption bands is associated with low water content. Therefore, the lower the R975/R900 ratio, the higher the water content of the Gorse companion plant. The two-way ANOVA showed a significant effect of the assemblage and water treatment. An interaction remained between these two factors (F = 5.709; p < 0.001). We thus compared the two water treatments (control and exclusion) for each assemblage. For the three species in the monospecific stand, the decrease in precipitation had no significant effects on the water contents compared to the controls (Figure 7). In the bispecific stands, Up improved the water-conserving abilities of Qc, while the opposite effect was observed for Ca. The Ro water contents were not affected by Up. However, when Qc was in a tri-specific stand with Ro or Ca, Up no longer improved the water content of Qc, nor did it even reduce it (QcRoUp). The water content of Ro in the trispecific stand was improved in the presence of Ca but reduced in the presence of Qc. In the tetraspecific stand, each of the three species had a better water capacity. High reflectance in water absorption bands is associated with low water content. Therefore, the lower the R975/R900 ratio, the higher the water content of the Gorse companion plant. The two-way ANOVA showed a significant effect of the assemblage and water treatment. An interaction remained between these two factors (F = 5.709; p < 0.001). We thus compared the two water treatments (control and exclusion) for each assemblage. For the three species in the monospecific stand, the decrease in precipitation had no significant effects on the water contents compared to the controls (Figure 7). In the bispecific stands, Up improved the water-conserving abilities of Qc, while the opposite effect was observed for Ca. The Ro water contents were not affected by Up. However, when Qc was in a tri-specific stand with Ro or Ca, Up no longer improved the water content of Qc, nor did it even reduce it (QcRoUp). The water content of Ro in the trispecific stand was improved in the presence of Ca but reduced in the presence of Qc. In the tetraspecific stand, each of the three species had a better water capacity.

Chlorophyll-a Fluorescence and Photosynthesis Index: NIRvP
We observed a significant interaction between the factor's assemblage and water treatment for the quantitative variable NIRvP: two-way-ANOVA (F = 11.123; p < 0.001). For each assemblage, the control and exclusion were thus compared (Figure 8). In the bispecific stand, Up enhanced photosynthesis in Ca as well as Ro under rain-exclusion devices. This could mean that Up increased the biomass of Ca and Ro under precipitationexclusion conditions. Furthermore, photosynthesis was not improved when Qc was associated with Up. Up therefore does not influence the biomass production in Qc under drought conditions. However, in the trispecific stand with Ca, Qc and Ro displayed an increase in photosynthetic capacity with rain exclusion.

Chlorophyll-a Fluorescence and Photosynthesis Index: NIRvP
We observed a significant interaction between the factor's assemblage and water treatment for the quantitative variable NIRvP: two-way-ANOVA (F = 11.123; p < 0.001). For each assemblage, the control and exclusion were thus compared (Figure 8). In the bispecific stand, Up enhanced photosynthesis in Ca as well as Ro under rain-exclusion devices. This could mean that Up increased the biomass of Ca and Ro under precipitation-exclusion conditions. Furthermore, photosynthesis was not improved when Qc was associated with Up. Up therefore does not influence the biomass production in Qc under drought conditions. However, in the trispecific stand with Ca, Qc and Ro displayed an increase in photosynthetic capacity with rain exclusion.

Chlorophyll-a Fluorescence and Photosynthesis Index: NIRvP
We observed a significant interaction between the factor's assemblage and water treatment for the quantitative variable NIRvP: two-way-ANOVA (F = 11.123; p < 0.001). For each assemblage, the control and exclusion were thus compared (Figure 8). In the bispecific stand, Up enhanced photosynthesis in Ca as well as Ro under rain-exclusion devices. This could mean that Up increased the biomass of Ca and Ro under precipitationexclusion conditions. Furthermore, photosynthesis was not improved when Qc was associated with Up. Up therefore does not influence the biomass production in Qc under drought conditions. However, in the trispecific stand with Ca, Qc and Ro displayed an increase in photosynthetic capacity with rain exclusion.

Nitrogen-Content Index and Other Vegetation Indices
The PLSDA analysis performed on all the selected indices (18) exhibited 4 hits with VIP scores >1 (Supplementary Figure S3): NPQI, NDVI(N), PRI and SIFb. The two-way ANOVA carried out with the nitrogen-content index, NDVI(N), showed a significant interaction between the assemblage and the water availability (F= 6.921, p > 0.001). For each Figure 8. Histograms of photosynthesis and fluorescence index (NIRvP) of assemblages in control (C) and rain-exclusion (E) conditions. Quercus coccifera (Qc), Cistus albidus (Ca) and Rosmarinus officinalis (Ro) and Ulex parviflorus (Up). The first two letters of each assemblage indicate the species measured. NS, not significant; * p < 0.05; ** p < 0.01; *** p < 0.001. n = 27.

Nitrogen-Content Index and Other Vegetation Indices
The PLSDA analysis performed on all the selected indices (18) exhibited 4 hits with VIP scores >1 (Supplementary Figure S3): NPQI, NDVI(N), PRI and SIFb. The two-way ANOVA carried out with the nitrogen-content index, NDVI(N), showed a significant interaction between the assemblage and the water availability (F= 6.921, p > 0.001). For each assemblage, the nitrogen levels were compared between the control and exclusion treatments. In the monospecific assemblages (Figure 9), the decrease in precipitation did not have a significant effect on the nitrogen content compared to the controls. However, in the bispecific assemblage, Up improved only the nitrogen content of Ca; in addition, it probably did so in the trispecific (CaQcUp) and tetraspecific (CaQcRoUp) assemblages. It should be noted that the nitrogen content was also increased under the drought conditions for QcCaUp and RoQcUp. assemblage, the nitrogen levels were compared between the control and exclusion treatments. In the monospecific assemblages (Figure 9), the decrease in precipitation did not have a significant effect on the nitrogen content compared to the controls. However, in the bispecific assemblage, Up improved only the nitrogen content of Ca; in addition, it probably did so in the trispecific (CaQcUp) and tetraspecific (CaQcRoUp) assemblages. It should be noted that the nitrogen content was also increased under the drought conditions for QcCaUp and RoQcUp. Figure 9. Histograms of nitrogen-content index, NDVI(N) of assemblages in control (C) and rainexclusion (E) conditions. Quercus coccifera (Qc), Cistus albidus (Ca) and Rosmarinus officinalis (Ro) and Ulex parviflorus (Up). The first two letters of each assemblage indicate the species measured. NS, not significant; * p < 0.05; ** p < 0.01; *** p < 0.001. n = 27.

Morphology and Anatomy of Ulex parviflorus
The SEM pictures of the Ulex parviflorus spines show that these were glabrous, with a high density of oval stomata ( Figure 10A,C). The stem shape was fluted, with abundant forked trichomes in pairs, but with a lower density of stomata compared to the spines ( Figure 10B

Morphology and Anatomy of Ulex parviflorus
The SEM pictures of the Ulex parviflorus spines show that these were glabrous, with a high density of oval stomata ( Figure 10A,C). The stem shape was fluted, with abundant forked trichomes in pairs, but with a lower density of stomata compared to the spines ( Figure 10B,D).
Plants 2023, 12, x FOR PEER REVIEW 1 assemblage, the nitrogen levels were compared between the control and exclusion ments. In the monospecific assemblages (Figure 9), the decrease in precipitation d have a significant effect on the nitrogen content compared to the controls. However, bispecific assemblage, Up improved only the nitrogen content of Ca; in addition, it ably did so in the trispecific (CaQcUp) and tetraspecific (CaQcRoUp) assemblag should be noted that the nitrogen content was also increased under the drought cond for QcCaUp and RoQcUp.

Morphology and Anatomy of Ulex parviflorus
The SEM pictures of the Ulex parviflorus spines show that these were glabrous a high density of oval stomata ( Figure 10A,C). The stem shape was fluted, with abu forked trichomes in pairs, but with a lower density of stomata compared to the s ( Figure 10B

Discussion
We provided a Gorse spectral signature, which may be distinguished from those of R. officinalis, Q. coccifera and C. albidus in terms of its architectural, structural and chemical organization. Our studies demonstrated how field spectrometry can discriminate between coexisting plant species. We showed that the spectral signature of Gorse differed from that of its companion species in specific spectral bands (Figure 3). The lower overall reflectance of the Gorse may be attributed to its spiny structure. However, the most discriminatory wavelength in the visible domain was exhibited by Rosmarinus and its assemblages ( Figure 5). The discrimination of Rosmarinus coexisting with Gorse and the two other species occurred mainly in green reflectance (559 nm). This clearly indicates that Rosmarinus' neighboring plants affect its chlorophyll-light-absorbing properties and, thus, photosynthesis. Our work agrees with other studies showing that hyperspectral measurements are applicable for taxonomic purposes [14]. This raises the question of the complex entanglement of the structural and/or molecular components involved in species distinction. We also showed that in natural conditions (control), the reflectance of individuals increases linearly and positively with the species richness. This is in accordance with a previous study carried out in the vicinity of the CLIMED site, which showed a linear increase in the plant productivity with species richness [32]. The occurrence of several kinds of species in the same area allows the exploration of different ecological niches, thus limiting interspecific competition. Indeed, the different root-distribution patterns of different species from the same community sharing the same ecological conditions is crucial for the understanding of how plants cope with their environmental constraints. The Q. coccifera roots and rhizomes were mainly located in the uppermost 15 to 35 cm of the soils [33], while the Gorse exhibited deep roots, with less lateral-root development in the upper layers [34]. In this study, both the kermes oak and the Gorse clearly explored different ecological niches. Therefore, the question is raised as to what would happen with a very high species richness. The linear profile observed would probably reach a plateau stage in accordance with the hypothesis of functional-trait redundancy.
In the context of climate aridification, the presence of Gorse also results in a modification of the spectral signature of its companion species, in this case, Cistus, rosemary and kermes oak ( Figure 6). However, after about ten years of excluding precipitation at the CLIMED site, we observed a total mortality of Gorse in monospecific stands under the exclusion devices. This clearly indicates a change in the landscape of the garrigue due to the ecological niches left vacant by the decline in the monospecific stands of Gorse, which raises many concerns. On one hand, the question of the species that would be most likely to occupy these new niches and, on the other hand, the mechanisms involved in the Gorse's decline. To address the latter question, the Gorse spines were analyzed through SEM images (Figure 10), which revealed a superficial arrangement of the stomata, with high density and devoid of protective appendages. This vulnerability to the increased atmospheric water demand in drought conditions would explain the great competition between Gorse individuals for the satisfaction of water needs. Indeed, Q. coccifera stomata are mainly located on the abaxial side, while R. officinalis and C. albidus have a high density of trichomes protecting their leaves from desiccation [17]. Gorse therefore persists in speciesdiverse environments, probably through niche diversity, which results in changes in the spectral response of its companion species. The most significant direct effect was observed with the Cistus (CaUp), whose reflectance in rain-exclusion conditions increased over the entire spectrum (Figure 6f; Supplementary Table S2). A higher reflectance means fewer internal-molecular-absorption features, as well as the involvement of structural aspects of leaf-like trichome abundance, wax thickness and texture [35]. This raises the question of the identification of the different molecules that may be characterized in such plantplant-microorganism interactions. Indeed, we identified wavelengths that discriminated the effect of drought and coexistence with Gorse (Supplementary Table S2). Each of these values could correspond to the radiation-absorption properties of specific compounds, leaf structures and canopy architectures. Recent studies have begun to identify the molecular intersections between biotic-and abiotic-stress responses through the concept of plant immunity [36]. Considering the overall richness of the communities investigated, we showed that drought induces a reduction in reflectance with the increase in the specific richness that follows a double reciprocal profile. Our data are in accordance with those of previous studies, which concluded that neighbors increase the profiles of local adaptation, with adaptation to competition occurring in benign habitats [37]. Further studies are required to characterize the role of niche diversity versus that of chemical metabolic changes.
In the specific case of the index linked to the water content of the tissues (ratio R975/R900), the effect of the Gorse varied according to the neighboring species. Gorse seems to exert a positive influence on kermes oak through its ability to conserve water in drought conditions and, therefore, has a facilitating role. This effect turns into antagonism in the case of Cistus, since it has a lower water content in the presence of Gorse in drought conditions. Because Gorse and kermes oak have different rooting systems, the observed facilitation effect is the result of niche diversity rather than the involvement of chemical mediators. However, the intense competition between Cistus and Gorse suggests that Gorse may act through the regulation of oak stomatal opening due either to volatile signal transducers or by the rhizosphere's chemical components. For instance, plants or rhizospheric microorganisms' exudation of abscisic acid is known to regulate stomatal conductance [38,39]. Considering all the tretraspecific assemblages, a decrease in the ratio of R975/R900 was noticed, confirming the same pattern as that obtained from the trispecific plant associations [17].
In the case of the role of Gorse as a Fabaceae for nitrogen supply in the ecosystem, a clear improvement was shown for Cistus in the bispecific, trispecific and tetraspecific stands. Although it is recognized as a releaser of allelopathic compounds, these do not affect the soil-nitrogen net rates of C. albidus [40]. Improvements were not observed in the bispecific assemblages of Q. coccifera and R. officinalis (Figure 9). For the rosemary and the kermes oak, the presence of nitrogen was increased in the trispecific stands, RoQcUp and QcCaUp, respectively. Our data raise the question of the nitrogen demand of wild species. Therefore, it is likely that Cistus' nitrogen metabolism is higher than that of the two other companion species. Nitrogen, as a macronutrient, is a part not only of the composition of proteins, but also of the structure of chlorophyll, which is a major pigment in photosynthesis [41]. Our work attempted to determine the possible variations in the photosynthetic capacity of species through the NIRvP index ( Figure 8). The U. parviflorus has a significant positive effect on the photosynthetic capacity of C. albidus. This is quite in accordance with the improvements in nitrogen content by the neighbor effect of Gorse, but also with other research reports [42,43]. A similar photosynthesis improvement was also obtained with Q. coccifera associated with C. albidus and Gorse. The R. officinalis NIRvP index increased when co-existing with the Gorse, but also in association with C. albidus. This result appears to have been dependent on the availability of water rather than nitrogen. In addition, the higher the number of species increases in the assemblages, the greater the photosynthetic capacity seemed, whether in natural conditions or in rain-exclusion conditions. Since NIRvP is positively correlated with gross primary productivity (GPP) and sun-induced fluorescence (SIF), our results indicate that NIRvP is a potential index to estimate the productivity of shrubland ecosystems in their specific diversity from space. Indeed, SIF is one of the major indicators for measuring the physiology of photosynthesis with remotesensing tools [44]. In addition, our work also identified other potential indicators of ecosystem structure and function, such as NPQI, PRI and SIFb.

Conclusions
Hyperspectral spectrometry has proven to be a promising nondestructive method for observing the neighbor effect in plant-plant interactions. Indeed, through the characterization of the spectral signature of Gorse, we followed the spatial dynamics of a shrubland following disturbances such as fires using remote-sensing devices. The spectroscopy also demonstrated that in rain-exclusion conditions, Gorse modifies the water content and the biomass of the three other companion species through photosynthesis, depending on their different assemblages. However, with the exception of Cistus, it does not modify nitrogen supply, as previously expected, even though it is a Fabaceae. Among the vegetation indices investigated, NPQI (normalized phaeophytinization index) appears as a robust canopy indicator of Gorse's neighbor effect. Spectroscopic measurements therefore make it possible to acquire valuable information on the physiologies of ecosystems independently of traditional field and laboratory surveys. This work also highlighted the vulnerability of Gorse to climate change, which urgently requires appropriate conservation measures. Therefore, a better understanding of biotic interactions seems necessary, particularly through the characterization of the metabolic networks underlying the observed spectral signatures.
Supplementary Materials: The following supporting information can be downloaded at: https://www. mdpi.com/article/10.3390/plants12040879/s1, Figure S1: Linear model of the integrated reflectance versus species richness in control condition. Figure S2: Correlation curve of the integrated reflectance versus species richness in rain exclusion condition. Table S1: Wavelengths (nm) that discriminate Ulex parviflorus from Quercus coccifera, Cistus albidus and Rosmarinus officinalis. ANOVA with p-value < 0.001 and post hoc analysis (Fischer's LSD test) from pairwise comparison in monospecific stands. Table S2: The most significant wavelengths between control (c) and rain exclusion treatment (e) for the different assemblages of Cistus albidus (Ca), Rosmarinus officinalis (Ro) and Quercus coccifera (Qc) in the presence of Ulex parviflorus (Up). ANOVA, p < 0.05. NS: non-significant.
Author Contributions: J.-P.M., acquired the funding, planned and designed the research, wrote the manuscript, and revised the manuscript; A.M. performed the reflectance measurements, wrote the manuscript, and revised the manuscript; M.F. performed the reflectance measurements, wrote the manuscript, and revised the manuscript. All authors have read and agreed to the published version of the manuscript.